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ABSTRACT 

Aims. The realism of current models of the penumbra is assessed by comparing their predictions with the plasma properties of 
penumbral filaments as retrieved from spectropolarimetric observations. 

Methods. The spectropolarimeter onboard Hinode allows us to distinguish for the first time the fine structure of the penumbra. 
Therefore, we can use one-component inversions to obtain the stratifications of plasma parameters in each pixel. The correlations 
between the plasma parameters and the continuum intensity are studied. 

Results. We find that, in the outer penumbra, the stronger flows and higher values of magnetic field inclination tend to be located in 
dark filaments. This finding does not seem to be compatible with the scenario of a field- free gappy penumbra. 
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1. Introduction 

There are two models of the penumbral fine structure under 
discussio n nowadays. The first is th e uncombed model pro- 
posed bylSolanki & Montavonl (|1993|) and further developed by 
iMartmez Pilletl (|2000|) to explain the broad-band circular polari- 
sation observed in sunspots. It envisages the penumbra as a col- 
lection of horizontal flux tu bes embedded in a mo r e verti cal and 
stronger background field. Schlichenm aier et aTl (|1998|) made 
simulations of the temporal evolution of such flux tubes. Initially 
placed at the sunspot magnetopause, they rise by magnetic buoy- 
ancy and quickly become horizontal in the photosphere. At the 
same time, a strong Evershed flow develops due to the gas pres- 
sure gradient that builds up during the rise. The tubes are hot- 
ter than the surrounding plasma and appear as bright penum- 
bral filaments which gradually cool down due to radiative losses, 
i.e. they become dark further from the bright penumbral grains 
where the tubes cross the r = 1 layer. Hereafter, this model will 
be called the rising flux tube model (RFT). 

Th e second model was suggested by [Spruit & Scharmen 1 
(2006). These authors argued that the RFT model cannot pro- 
duce sufficient heating of the penumbra. As an alternative, they 
proposed a field-free gap model (FFG) which consists of field- 
free material protruding into the penumbral magnetic field from 
below. With such a configuration, the bright filaments would be 
heated all along their lengths and not only at the point of emer- 
gence as in the case of the RFT model. While the FFG model 
does not offer any explanation for the Evershed flow, there is 
only one place where the flow can reside, namely on top of the 
field-free gaps, since this is the sole region in the model where 
horizontal fields exist. 

Both models have problems fro m a theoreti c al poin t of view 
which need to be clarified (see Bell ot Rubia |2007[ and ref- 
erences therein). However, as illustrated in Fig. [T] RFTs and 
FFGs should produce similar stratifications of the plasma pa- 



rameters if they are placed near r = 1 (marked by the gray 
line). The two models predict weaker and more inclined fields 
low in the atmosphere, along with increased line-of - sight (LOS) 
velocities. Using Hinode spectropolarimetric data, I Jurcak et aD 
(120071) found this configuration in the bright filaments of the in- 
ner penumbra. 

In this Letter we concentrate on the mid and outer penum- 
bra, where differences between the models could be more pro- 
nounced. According to the FFG model, horizontal fields and 
increased LOS velocities should be associated with bright ar- 
eas everywhere in the penumbra. By contrast, the RFT model 
expects such characteristics in the dark filaments of the outer 
penumbra. 

2. Observations and data analysis 

The data analysed here were obt ained using the spec- 
tropolarimeter (SP; iTarbell et all [2007) of the Solar Optical 




Fig. 1. Sketch of a RFT (left) and a FFG (right) surrounded by 
the background magnetic field of the penumbra. The gray line 
represents the possible position of the optical depth unity layer. 
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Fig. 2. Maps of continuum intensity at 630.3 nm reconstructed 
from the SP raster scans. North is up and West to the right. The 
maps cover the limb side part of the penumbra in evolving AR 
10930. The arrows point to the disc centre. The black and white 
bands indicate the areas selected for analysis. 



the light of Hinode spectropolarimetric observations 
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Fig. 3. Scatter plots of LOS velocity (left) and inclination (right) 
fluctuations vs. continuum intensity fluctuations. The black, red, 
and blue symbols represent points from the M, B, and O bands, 
respectively. 



Telescope (iTsuneta et all I2007I) onboard the Hinode satellite 
(Kosugietal ., 2007). This instrument measures the Stokes pro- 
files of the two iron lines at 630.15 nm and 630.25 nm. 

Normal SP scans of AR 10930, providing a spatial resolution 
of 0732 and a noise level of 10" 3 / c , were taken on December 
13, 14, and 15, 2006, when the spot was located at heliocen- 
tric angles of 30°, 41°, and 54°, respectively. The three panels 
of Fig. [2] show the limb-side penumbra of AR 10930 as recon- 
structed from the continuum intensities observed redward of the 
630.25 nm line. The Stokes V profiles emerging from the penum- 
bra exhibit significant area asymmetries due to the existence of 
gradients of the plasma parameters in the line-forming region. 
We can clearly see the bright and dark filaments in the contin- 
uum maps. Their signals are not mixed as in the case of ground- 
based observations, implying that the penumbral fine structure is 
distinguished in spectropolarimetric measurements for the first 
time. Higher spatial resolution observations of the penumbra ex- 
ist (e.g., Scharmer et al. 2002; Rouppe van der Voort et al. 2004; 
Rimmele 2004; Langhans et al. 2005; Bellot Rubio et al. 2005; 
Sanchez Almeida et al. 2007), but they correspond to filtergrams, 
Dopplergrams, magnetograms, and spectroscopic measurements 
from which the magnetic field vector cannot be determined. 

Each day we selected three 5 -pixel- wide azimuthal cuts. As 
illustrated in Fig. [2 they sample the mid and outer parts of the 
penumbra. The bands were chosen so as to minimise any trend 
of the continuum intensity along them, except for small-scale 
fluctuations caused by the alternating bright and dark filaments. 
In total, they contain 5862 pixels. 

Th e data reduction is described in detail in Jurca k et aD 
d2007h . The zero point of the velocity scale is taken to be the line 
core position of the mean quiet Sun intensity profile, computed 
separately for each slit. The observ ed Stokes spectra ha ve been 
inverted using the SIRGAUS co de dBellot RubioL |2003 " |). which 
is a modified version of SIR dRuiz Cobo & del Toro IniestaL 
119921) . This code presumes the existence of a Gaussian perturba- 
tion (GP) in the stratifications of plasma parameters somewhere 
in the line-forming region. Given the high angular resolution 
of the Hinode measurements, we only consider one-component 
model atmospheres. 

The inversion code looks for the best solution in a space of 
13 free parameters. Six of them define the physical conditions of 
the unperturbed atmosphere: we use two for the temperature (T), 
and one for the field strength (B), inclination (yix>s), azimuth 
(0los), and LOS velocity (vlos)- The width and position of the 
GP is the same for all plasma parameters, which adds two more 
free parameters. The rest of parameters are the amplitudes of 
the perturbation in the different physical quantities. As initial 



values we adopt AT = +800 K, AB = -500 G, Ay L os = -30°, 
A0los = -5°, and Avlos = +3 km s" 1 . These initial parameters 
imply weaker and more inclined fields associated with stronger 
Evershed flows, so they can model both RFTs and FFGs. The 
magnetic field inclination and azimuth retrieved by the inversion 
code are line-of- sight values. We transformed these variables to 
the local reference frame. In what follows, the field inclination 
(y) represents the angle between the magnetic field vector and 
the inward normal to the solar surface. 



3. Choice of plasma parameters 

As explained by Jurca k et all (l2007|) . we perform four inversions 
per pixel with different starting heights of the GP, and select the 
one that shows the smallest value of the merit function (sum 
of squared differences between observed and synthetic profiles). 
It sometimes happens that the Stokes spectra emerging from a 
given pixel can be fitted equally well by GPs located either in 
deep layers or in the mid photosphere. In these cases it is dif- 
ficult to decide which solution is closer to reality on the ba- 
sis of the merit functions alone, as they are usually of compa- 
rable value. Although the resulting stratifications are different, 
they share some similarities: variations in the position of the 
GP mainly change the amplitude and sign of the magnetic field 
strength perturbation, but the sign and amplitude of the GP for 
the LOS velocity and the field inclination are almost independent 
of the height of the Gaussian. 

We take advantage of this fact by adopting the part of the so- 
lution which is similar in all cases, i.e., the average values of the 
plasma parameters in the line-forming region. Specifically, we 
compute for each pixel the mean value of the LOS velocity, field 
inclination, and field strength in the range from logr = -0.5 to 
logr = -2.5. This is the region where the Stokes profiles of the 
Fe i 630.2 nm lines are more sensitive to the physical conditions 
of the atmosphere (Cabrera Solana et al. 2005). 

From the obtained values we subtract the locally averaged 
plasma parameters, as defined by a running mean over 100 pix- 
els. The resulting fluctuations of continuum intensity, LOS ve- 
locity, field inclination, and field strength are used to test the 
predictions of the RFT and FFG models. 

4. Results 

Figure [3] shows scatter plots of LOS velocity and magnetic 
field inclination fluctuations vs. continuum intensity fluctua- 
tions. There is a clear tendency for weaker flows and more ver- 
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Fig. 4. Continuum intensity (black), LOS velocity (blue), and magnetic field inclination (red) fluctuations along the selected paths. 
The plots show, from top to bottom, parts of the Ml, Bl, M3, and B3 bands. 



tical fields in brighter penumbral areas. However, this type of 
plot does not reflect how the plasma parameters are associated 
with penumbral filaments. Therefore, in Fig.|4]we give examples 
of continuum intensity fluctuations (black lines) and the corre- 
sponding changes of LOS velocity (blue lines) and inclination 
(red lines) along one-pixel-wide azimuthal cuts extracted from 
the selected bands. 

Table Q] summarises the correlation coefficients obtained in 
the nine bands, with errors representing 95% confidence inter- 
vals. The values shown in the second and third columns confirm 
the visual impression from Figs. [3] and 01 that, in the mid and 
outer penumbra, the stronger flows and the more inclined fields 
usually occur in dark filaments. The average correlation coeffi- 
cient between Avlos and A/ c amounts to -0.40, which is quite 
high compared to previous analyses based on spectropolari - 
metric data (e.g. JWestendorp Plaza etall 1200 ll: ISolankiL l2QQ3h . 
Enhanced LOS velocities in dark filaments were also found from 
high-resolut ion filtergrams with comparable or even better corre - 
lations (e.g. JSchlichenmaier et all 120051: llchimoto et alll2007l) . 

The third column of Table Q] shows the correlation coeffi- 
cients between Ay and A/ c . In this case the average correlation 
is -0.40 (-0.46 if we do not take into account bands Bl and 
01, which are significantly worse). Thus, there is a tendency for 
the more horizontal fields to be associated with dark penumbral 
areas. Previous analyses reported either no signific ant correla- 
tion between these parameters i n the outer penumbra dTitle et all 

or larger field inclinations in 
dark filaments (e.g., [ Rimmele, 1 19951: IWestendorp Plaza et all 
1200 ll: lLanghans et all l2005h . Usually, the correlation coeffi- 
cients quoted in those works are smaller (in absolute value) than 
0.2. Only when the correlation is analysed at small spatial scales 
do the coefficients become similar to (or higher than) the ones 
reported here. 



Table 1. Correlation coefficients between fluctuations of LOS 
velocity and continuum intensity (second column), field inclina- 
tion and continuum intensity (third column), field inclination and 
LOS velocity (fourth column), and field inclination and strength 
(fifth column), along the nine bands marked in Fig. 2. Errors 
represent 95% confidence intervals. 





Avlos 


Ay 


Ay 


Ay 


Cut 


vs. A/ c 


vs. A/ c 


vs. Avlos 


vs. AB 


Ml 


-0.56 + 0.06 


-0.55 + 0.06 


0.82 + 0.03 


-0.78 + 0.04 


M2 


-0.23 + 0.08 


-0.44 + 0.07 


0.78 + 0.03 


-0.75 + 0.04 


M3 


-0.44 + 0.07 


-0.40 + 0.07 


0.78 + 0.03 


-0.69 + 0.04 


Bl 


-0.34 + 0.08 


-0.15 + 0.09 


0.81+0.03 


-0.72 + 0.04 


B2 


-0.36 + 0.07 


-0.53 + 0.05 


0.67 + 0.04 


-0.71+0.04 


B3 


-0.41 + 0.07 


-0.43 + 0.06 


0.70 + 0.04 


-0.46 + 0.06 


01 


-0.30 + 0.07 


-0.18 + 0.07 


0.86 + 0.02 


-0.30 + 0.07 


02 


-0.42 + 0.06 


-0.44 + 0.06 


0.77 + 0.04 


-0.33 + 0.07 


03 


-0.55 + 0.05 


-0.45 + 0.06 


0.60 + 0.04 


-0.30 + 0.06 



The two last columns of Table \T\ quantify the relation be- 
tween the changes of magnetic field inclination and the LOS 
velocity and field strength. It is generally accepted that the 
Evershed flow resides in areas of horizontal fields. With cor- 
relation coefficients as high as 0.7-0.8, the fourth column of 
Table [T] just confirms this fact. The fifth shows that the more 
inclined fields are weaker, although the correlation coefficients 
decrease significantly from the mid to the outer penumbra. This 
implies a smaller difference in magnetic field strength between 
areas with horizontal fields (or flow channels) and the surround- 
ing atmosphere in the outer penumbra, which is not surprising 
as both the RFT and FFG models predict an equalisation of the 
field strength with increasing distance from the umbra. A smaller 
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field strength difference between inclined flow channels and sur- 
rounding plasma has also been inferred from two-component 
Stoke s inversions (e.g.. iBellot Rubio et all 120041; iBorrero et all 
2006) and analyses of the net circular p olarisation observed i n 
the penumbra at high angular resolution dTritschler et alll2007l) . 

The plasma parameter averages made over the line-forming 
region do not tell us much about the actual properties of the 
flow channels represented by the GP. Therefore, we also stud- 
ied the general behaviour of the peak values of the GPs of 
LOS velocity and magnetic field inclination. We found that 
the flow becomes faster with increasing distance from the um- 
bra (on average it is 2 km s" 1 faster in the O cuts than in 
the M cuts). This increase confirms previous reports of en- 
hance d flow ve l ocity with increasing distance from the umbra 
(e.g. iRimmelei 119951: IWestendorp Plaza et all l200lt ISolankiL 
12003 : Bellot Rub io et all 12004 IBorrero et all l2006h . and is 
also predicted by simulati ons of moving penumbral tubes 
(ISchlichenmaier et allfl998l) . 

The maximum values of the LOS velocity are around 
10 km s" 1 near the symmetry line in the 02, B3, and 03 cuts. 
Even if such velocities are overestimated because they represent 
the peak amplitudes of the GP, the absolute flow velocity must 
be very high due to projection effects. These velocities are com- 
patible with the supersonic flows predicted by the simulations of 
Schli chenmaier et all (Il998l) in the outer penumbra. 

The average magnetic field inclinations in the flow channels 
are 85°, 100°, and 105° in the M, B, and O cuts, respectively. 
This means that the flow channels are inclined upward in the mid 
penumbra and point down in the outer penumbra (the polarity 
of the spot is negative). Magnetic fields returning to the solar 
surface in the outer pe numbra we re already reported by, e .g., 
Westendorp Plaza et al . (2001) and Bellot Rubio_etaL| <2004). 

5. Discussion and conclusions 

The most important result of our analysis is that we do not find 
positive correlation coefficients between I c and y or vlos- This 
indicates that bright penumbral filaments show the more verti- 
cal fields and weaker flows, in contradiction with the predictions 
of the FFG model. The horizontal fields and stronger flows tend 
to be associated with dark filaments, as has been reported pre- 
viously from lower resolution data. The fact that the correlation 
coefficients between I c and y or between I c and vlos are not - 1 
can be easily explained by the mixing of flux tubes and back- 
ground magnetic fields. In the mid and outer penumbra, dark 
areas may represent either the darkened tails of RFTs or back- 
ground fields. The latter are characterised by more vertical fields 
and weaker flows. Such an opposite behaviour decreases the 
correlation coefficients, which nevertheless remain negative be- 
cause of the overall dominance of dark RFTs in the outer parts of 
the penumbra. This argument does not influence the correlation 
coefficients between y and vlos, which are indeed much higher 

tha n those involving Z c . 

ISchlichenmaier et al J (120051) and llchimoto et al.l (120071) came 
to a similar conclusion (regarding the velocities) from sim- 
pler interpretations of filtergrams and spectropolarimetric data, 
respectively. They also reported that the horizontal fields and 
the flows are assoc iated with bright filaments in the innermos t 
penumbra (see also iHirzberg er et al. i 120051: Uurcak et all [2007), 
where the dark areas most likely represent background magnetic 
fields. Thus, the correlation coefficients between continuum in- 
tensity and LOS velocity/inclination change sign from the inner 
to the outer penumbra, due to the fact that RFTs predominantly 
show up as bright structures in the inner penumbra but darken 



toward the outer edge of the spot. This sign reversal cannot be 
explained by the FFG model. The existence of magnetic field 
lines pointing to the solar interior in the outer penumbra, also 
confirmed by o ur observations, repr esents another problem for 

the FFG model dBellot RubiqU2007^ 

The RFT simulations of ISchlichenmaier et al.l ([T998) seem 
to be able to explain all the results of our analysis, including 
the sign reversal of the correlation coefficients and the larger 
mean value of the LOS velocity in the outer penumbra as com- 
pared with the mid penumbra. Overall, this suggests that the RFT 
model provides a better description of the penumbra than the 
FFG model. 

Acknowledgements. This work has been enabled thanks to the funding provided 
by the Japan Society for the Promotion of Science. Hinode is a Japanese mission 
developed and launched by ISAS/JAXA, with NAOJ as domestic partner and 
NASA and STFC (UK) as international partners. It is operated by these agencies 
in cooperation with ESA and NSC (Norway). The computations were carried out 
at the NAOJ Hinode Science Center, which is supported by the Grant-in- Aid for 
Creative Scientic Research The Basic Study of Space Weather Prediction from 
MEXT, Japan (Head Investigator: K. Shibata), generous donations from Sun 
Microsystems, and NAOJ internal funding. Financial support from the Spanish 
Ministerio de Education y Ciencia through project ESP2006-13030-C06-02 is 
gratefully acknowledged. 



References 

Bellot Rubio, L. R. 2003, ASP Conf. Series, 307, 301 
Bellot Rubio, L. R., Balthasar, H., & Collados, M. 2004, A&A, 427, 319 
Bellot Rubio, L. R., Langhans, K., & Schlichenmaier, R. 2005, A&A, 443, L7 
Bellot Rubio, L.R. 2007, in: Highlights of Spanish Astrophysics IV, eds. F. 

Figueras et al. (Dordretch: Springer), 271 
Borrero, J. M., Solanki, S. K., Lagg, A., S ocas -Navarro, H., & Lites, B. 2006, 

A&A, 450, 383 

Cabrera Solana, D., Bellot Rubio, L.R., & del Toro Iniesta, J.C. 2005, A&A, 
439, 687 

Hirzberger, J., Stangl, S., Gersin, K., et al. 2005, A&A, 442, 1079 
Ichimoto, K., Shine, R.A., Lites, B.W., et al. 2007, PASJ, 59, 593 
Jurcak, J., Bellot Rubio, L.R., Ichimoto, K., et al. 2007, PASJ, 59, 601 
Kosugi, T., Matsuzaki, K., Sakao, T., et al. 2007, Sol. Phys., 243, 3 
Langhans, K., Scharmer, G. B., Kiselman, D., Lofdahl, M. G., & Berger, T. E. 

2005, A&A, 436, 1087 
Martinez Pillet, V. 2000, A&A, 361, 734 
Rimmele, T. R. 1995, A&A, 298, 260 
Rimmele, T. R. 2004, ApJ, 604, 906 

Rouppe van der Voort, L. H. M., Lofdahl, M. G., Kiselman, D., & Scharmer, 

G. B. 2004, A&A, 414, 717 
Ruiz Cobo, B. & del Toro Iniesta, J. C. 1992, ApJ, 398, 375 
Sanchez Almeida, J., Marquez, L, Bonet, J. A., & Dommguez Cerdena, I. 2007, 

ApJ, 658, 1357 

Scharmer, G. B., Gudiksen, B. V., Kiselman, D., Lofdahl, M. G., & Rouppe van 

der Voort, L. H. M. 2002, Nature, 420, 151 
Schlichenmaier, R., Jahn, K., & Schmidt, H. U. 1998, A&A, 337, 897 
Schlichenmaier, R., Bellot Rubio, L. R., & Tritschler, A. 2005, AN, 326, 301 
Solanki, S. K. & Montavon, C. A. P. 1993, A&A, 275, 283 
Solanki, S. K. 2003, A&A Rev., 11, 153 
Spruit, H. C. & Scharmer, G. B. 2006, A&A, 447, 343 
Stanchfield, D. C. H., II, Thomas, J. H., & Lites, B. W. 1997, ApJ, 477, 485 
Tarbell, T., et al. 2007, Sol. Phys., in preparation 

Title, A. M., Frank, Z. A., Shine, R. A., Tarbell, T. D., Topka, K. P., Scharmer, 

G., & Schmidt, W. 1993, ApJ, 403, 780 
Tsuneta, S., Ichimoto, K., Katsukawa, Y., et al. 2007, Sol. Phys., submitted 
Tritschler, A., Muell er, D. A. N„ S chlichenmaier, R., & Hagenaar, H. J. 2007, 

ApJ, in press, larXiv:07 10.45451 
Westendorp Plaza, C, del Toro Iniesta, J. C, Ruiz Cobo, B., & Martinez Pillet, 

V. 2001, ApJ, 547, 1148 



